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Abstract 

The magnitude of the strong interaction is characterized by a s , the coupling parameter in Quan- 
tum Chromodynamics (QCD), a parameter with an unexplained value in the Standard Model. In 
this paper, a candidate explanation for a s is derived from (1) the lifetime of quark-antiquark pairs 
in vacuum fluctuations given by the Uncertainty Principle, (2) the variation of a s as a function 
of energy in QCD, and (3) classical relativistic dynamics of the quarks and antiquarks. A semi- 
classical model for heavy quark-antiquark vacuum fluctuations is described herein, based on (2) 
and (3). The model in this paper predicts the measured value of a s (M z o) to be 0.121, which is in 
agreement with recent measurements within statistical uncertainties. 
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I. INTRODUCTION 



As a result of the Uncertainty Principle, vacuum fluctuations occur consisting of a par- 
ticle and its antiparticle created by the vacuum and annihilated in a short lifetime At. A 
fluctuation consisting of a quark and its antiquark interacting via a gluon can occur in per- 
turbation theory in QCD j^j], J3I), the quantum field theory that successfully describes 
quark interactions with great accuracy. 

For educational purposes, the author developed a classical model for the dynamics of such 
a quark- antiquark pair - a model that was intended to contrast the quantum mechanical 
prediction of the pair lifetime with the classical prediction. A surprising result was found: 
the pair lifetime in the classical model agrees to a good approximation with the lifetime in 
quantum mechanics. 

The model is described completely in this paper. The model becomes semiclassical in the 
usual sense in a natural way. Because the pair lifetime from the model agrees so precisely with 
the quantum mechanical lifetime, particularly in the case of bottom quarks, the expressions 
for the lifetime in the two theories enable us to solve for the QCD strength parameter a s 
solely from these theoretical considerations. 

In the Standard Model, QCD is a physical asymptotically free field theory for any appro- 
priate value of the input parameter A, which is set via experimental measurement. Within 
the QCD theory alone it would be impossible to establish a s without recourse to a measure- 
ment. But the Uncertainty Principle and semiclassical QCD are not one and the same, and 
because we find in this work that the Uncertainty Principle and semiclassical QCD come 
to agreement on the vacuum fluctuation lifetime, we can logically derive what the physical 
measurement must be. It is shown below that this enables us to ground a s on the value of 

h. 



II. THE MODEL BASIS ENERGY CONSERVATION 

If a quark and its antiquark are positioned at rest in their center of mass reference frame 
and are released, then in general their mutual attraction will draw them to a collision at the 
origin where they will annihilate. Photons or other particles would result, given sufficient 
energy. However, the quark and antiquark experience a potential energy U(R) [4| as a 
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FIG. 1: Feynman diagram of a virtual quark-antiquark pair (VQAP). Time increases from bottom 
to top of figure. The curling segment represents gluon interaction. 

function of their separation distance that could reduce the mass-energy of the system, if the 
separation R between particles is small enough. The energy conservation relation 

e = 2-im q c 2 + U(R) =0 (1) 

is possible to satisfy, where the first term is the total relativistic energy of the particles, 
kinetic plus rest mass. (R = 2r with r the radius of either particle from the center of mass.) 
In classical physics, the collision would not yield any energy and so no photons or particles 
could be emitted. 

If the time-reversed ballistic trajectory occurred, with the vacuum spontaneously creating 
a quark-antiquark pair obeying Eq. (GQ), then the particles only could move apart in one- 
dimensional motion to reach turning points separated by R m ax- Continuing this trajectory 
so that the particles fall from the turning points back to the origin, they would disappear 
back into the vacuum, like a virtual quark-antiquark pair (VQAP; see Fig.dJfor the Feynman 
diagram). This is similar to virtual electron-positron pairs, as discussed by Greiner (p. 3 of 



ref. j^]) and Sakurai (p. 139 of ref. jf]]). 



Quantum theory implies that this two-particle system of quarks would obey the time- 
energy uncertainty relationship for the energy fluctuation Ae (p. 139 of ref. js]) 

AeAt = = 5.273 x l(T 28 erg s. (2) 

The energy fluctation is Ae = 2m q c 2 , since the mass-energy of each quark contributes m q c 2 . 
The quantum-mechanical lifetime of the fluctuation is At. Since h is the quantum of action, 
Eq. (J2J) establishes an action integral A that characterizes a VQAP that has Ae = 2m q c 2 . 

The first purpose of this paper is to present classical computations of At and the action 
integral for the trajectory described above, which turn out to give results that satisfy Eq. (j2J) 



remarkably well, provided that the QCD interaction between the particles is well-described 
by the potential energy function. 

The second purpose of the paper follows from the fact that QCD cannot specify the value 
of a s at arbitrary energy or 4-momentum scale Q without an established measurement of a s 
at some particular energy \i [7|; but once the renormalized coupling a s (/i 2 ) is measured, then 
QCD precisely gives the variation of a s as a function of energy (the "running" coupling). 
The present paper offers a theory based on the action integral that establishes the value of 
a s at the energy scale of twice the bottom quark mass-energy to good approximation. This 
enables one to use the QCD running coupling a s (Q 2 ) to determine the coupling strength in 
general in the usual way to good approximation. 



III. QCD POTENTIAL ENERGY FUNCTION 

As discussed in detail by Lucha et al. (especially pp. 161-162 of ref. a potential 
energy function serves to describe the bound states of heavy quarks (charm, bottom, and 
top). For light quarks the QCD interaction is not satisfactorily described by a potential 
energy function and will not be attempted here. Here we apply the standard potential 
energy treatment to the heavy charm, bottom and top quarks. 

We use the standard Cornell potential sj], js] 

V(R) = -p^ + aR, (3) 

where a s is the dimensionless QCD strong coupling strength and a ~ 0.25 GeV 2 . The second 
term, aR, is only significant for R > 10~ 13 cm. We will not need to consider the aR term, 
since the first term with the Coulomb-like dependence turns out to strongly dominate the 
potential because R ma x "C 1CT 13 cm for VQAPs. 

The model VQAP is a form of bound state. The standard way to account for the variation 
of a s (Q 2 ) in a quark bound state is to let a s depend on the quark masses and use Q 2 = 
(mi + 1712) 2 , with the rrii the quark masses (see p. 129 of ref. 10]). To model a VQAP we 
may then compute a s to leading order (Eq. (6) of Ref. [?}|). So the QCD coupling is then 
given by 
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where /3q is defined by 



_ 33 - 2n f 



rif is the number of quark flavors with masses much less than mi + m 2 , and A is the QCD 
scale energy, 



2 



A 2 = E (Q) 



e 



A is approximately 0.093 GeV, assuming that \i = M^o = 91.2 GeV (the mass-energy of 
the Zq particle), a s (M z o) = 0.119 ± 0.002. This is a typical value of A for a one-loop 
approximation ([7] p. R31). 

For the running coupling near the charm quark mass rif = 4, near the bottom quark mass 
rif = 5, and near the top quark mass rif = 6. We now can use ni\ = ni2 = m q for each quark 
mass m q in the expression for Q 2 and compute a s from Eq. (jlj) for a running coupling. Then 
with a charm quark current mass m q = m c of 1.27 GeV/c 2 we have a s (2m c c 2 ) = 0.228; with 
a bottom quark current mass m q = vn,^ of 4.2 GeV/c 2 we have a s (2rribC 2 ) = 0.167; with a 
top quark current mass m q = m t of 173.1 GeV/c 2 we have a s (2m t c 2 ) = 0.109 



IV. CLASSICAL BALLISTIC TRAJECTORY AND LIFETIME OF VQAP 

Let us calculate the trajectory lifetime t vq for a VQAP. If one solves the energy equation 
([T|) for the dynamics using V(R), the nonrelativistic potential energy, then the particle 
velocities nevertheless exhibit relativistic motion, approaching c asymptotically at the origin 
r = 0. Thus it is necessary to correct the potential energy for relativistic effects. Jackson 
([l^. p. 553) demonstrates how this is done by using the relativity factor 7 defined in 
Eq. (TlOl) below. For this trajectory of linear motion, the transformation R — > 7/? in the 
expression for V(R) performs the appropriate modification of the potential energy function 
(since we are considering the center-of-mass reference frame). The potential energy function 
in Eq. ([3]) becomes 

U(R) = -ii He. (7) 

With this U(R) we can solve Eq. (TTJ) for R ma x at the turning point (where 7=1): 

4 a 

2m q c 2 = - U(R max ) = -— —he (8) 

R = (9) 
3 m q c 



For the charm quark current mass m c of 1.27 GeV/c 2 = 2.26 x 10 24 g, we find the charm 
VQAP has an R max = 2.35 x 10" 15 cm. 

Checking the terms in Eq. fl3]) for V(R max ) shows that the first term is about -100 times 
the second term. This confirms that the quarks are so deep in the potential well that the 
aR term of the Cornell potential can be neglected in solving the problem. 

Let us define the time from appearance of a quark at the origin r = to the time that 
the quark stops at the turning point r = \R m ax as \t vq . The t vq is the classical equivalent 
of the quantum-mechanical At that we seek. We note that 

and solve for dt, which we shall integrate. We rewrite the energy equation flTJ with ( = 

Rj Rmax aS 

f^**-^. (ID 

The time for the particle to fall from r = R max /2 back to r = is also ht vq , so we have 

rR max /2 dr R max r 1 d( R max y^T(l) _ 2R max _ 4 a s h 
vq ~ Jo Cv^C" c Jo V^C~ c r(|) c ~3m 9 c 2 1 j 

The integral is given in ref. ( 13j, p. 974). The value of t vq is the total time for either quark 
to travel from r = to its turning point and back to r = 0. 

For the charm quark, with m q = m c fa 2.26 x 10 -24 g, we find the trajectory lifetime of the 
VQAP to be t vq f« 1.57 x 10 -25 s. In comparison, the standard lifetime of the charm VQAP, 
given by the Uncertainty Principle expressed in Eq. (jSJ), is At = h / (4m c c 2 ) ~ 1.29 x 10~ 25 
s. So t vq from the classical computation is approximately 22% larger than At. This is 
remarkably close agreement of the classical lifetime with the quantum mechanical lifetime. 

For the bottom quark, with m q = rrib ~ 7.48 x 10 -24 g, we find the trajectory lifetime of the 
VQAP to be t vq ~ 3.47 x 10 -26 s. In comparison, the standard lifetime of the bottom VQAP, 
given by the Uncertainty Principle expressed in Eq. (121), is At = h/(4mf,c 2 ) ~ 3.90 x 10~ 26 
s. So t vq from the classical computation is approximately 11% smaller than At. This also is 
remarkably close agreement of the classical lifetime with the quantum mechanical lifetime. 

For the top quark, with m q = m t ~ 3.08 x 10~ 22 g, we find the trajectory lifetime of the 
VQAP to be t vq fa 5.51 x 10 -28 s. In comparison, the standard lifetime of the top VQAP, 
given by the Uncertainty Principle expressed in Eq. (j2J), is At = h/(4m t c 2 ) fa 9.48 x 10 -28 



6 



s. So t vq from the classical computation is approximately 42% smaller than At. This also is 
remarkably close to the quantum mechanical lifetime. 

V. ACTION INTEGRAL FOR THE TRAJECTORY 

A key step in quantizing a classical model to make it a semiclassical model of a quantum 
system is computation of the action integral. In the present model, that is done as follows. 
The expression for the integral of action associated with a potential function U acting on a 
particle, in the relativistic case, is given by Lanczos (|14j. p. 321): 

A=-P (13) 

Jti c 

where ds = cdt/'y. Considering the integrated action of the potential energy field in a 
VQAP, we compute the field action integrated over t vq : 

rtvq/2 / 4a s fic\dt 8a s fi /" R ™^/ 2 dr 
4a s h r 1 d( 4a s h r 1 d( 8a 



3 Jo C-KVT^C 3 Jo VT^C 3 H (15) 
For the charm VQAP, a s (2m c c 2 ) = 0.228 and therefore A = 0.61ft. This action integral is 
only 22% larger than the exact VQAP quantum fluctuation action in Eq. (J2]), \fi- 

In the case of the bottom quark, the action integral for the model of the VQAP is found 
by substituting a s (2mbC 2 ) = 0.167 into Eq. ( 115]) . and we find A = 0.45 This is 10% lower 
than the quantum mechanical action for a VQAP. 

In the case of the top quark, the action integral for the model of the VQAP is found by 
substituting a s (2m t c 2 ) = 0.109 into Eq. (JTJ, and we find A = 0.29 ft. This is 42% lower 
than the quantum mechanical action for a VQAP. 

VI. SPIN-SPIN INTERACTION EFFECTS 



As noted by Lichtenberg (p. 133 of ref. 10]) the spin-spin interaction between quarks 



is the source of electromagnetic mass splittings among hadron isospin multiplets. Here 
it will be shown that spin-spin interaction between the quark and antiquark in a VQAP 
is significant for the charm and top quark cases. Spin-spin interaction would modify the 
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dynamics of the semiclassical model and lead to a different action integral and trajectory 
lifetime, so its influence needs to be quantified. 

The creation of a quark-antiquark pair from the vacuum should conserve angular mo- 
mentum, so we consider the case in which the particles have spin parallel to the motion axis 
but pointing in opposite directions. The potential energy function in Eq. (J7J) becomes ( 12 1 

(16) 



where 



gQ q h 



* = 2^2 (17) 
is the magnetic moment of a quark with electric charge Q q and g is the Lande factor, which 
will be taken as equal to 2 for present purposes. 
The energy equation Eq. (JTJ) in this case becomes 

£ = W-~Sc-^ = (18) 

Again we can derive the separation of the particles at the turning point R' max by letting 
7 — > 1. It is convenient to use ( = R/R max as a dimensionless variable again, and with that 
substitution, the energy equation becomes 

270 2 

Cl,-C-A = with A = ^-. (19) 

The parameter A is a measure of the magnitude of the spin-spin interaction relative to the 
QCD interaction; A — > recovers the previous case in Eq. ([I]). 

Solutions of this new energy equation (1TB]) are deferred for a later paper, but the values 
of A for each of the heavy quarks are illuminating. The electric charges of the charm and 
top quarks are both |e, but the charge of the bottom quark is only — |e . Consequently the 
parameter A is 0.23 for the charm quark and 2.1 for the top quark. For the bottom quark, 
A = 0.15, which indicates that the importance of spin-spin effects in the dynamics of the 
model is minimal for the bottom quark. 

The smallest value of A is associated with the quark that exhibits the best agreement 
between At and t vq , the bottom quark. This suggests that the bottom quark case is best 
for using the model to characterize the QCD interaction without perturbations from the 
spin-spin effects. 
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In the above cases, the model's representation of the bottom VQAP inherently is approx- 
imately quantized - a remarkable agreement between a quantum-mechanical characteristic 
of a dynamical system and the classical description of it. In comparison, semiclassical mod- 
els for mesons, which achieve excellent agreement with measurements of meson masses [l^] , 
[3| need to be formulated with additional quantization conditions that introduce the factor 
h. We have not imposed any quantization conditions upon the trajectory in this dynamical 
model. The bottom quark model herein achieves approximate quantization at a s (2rribC 2 ) 
based upon only the measured value of a s (M z o), the QCD theoretical energy dependence 
of a s (Q 2 ), and relativistic dynamical theory (Eqs. ([rj and (J7D). 

Preliminary work to investigate the discrepancy between t vq and At in the cases of the 
charm quark and top quark has been performed as the author will show elsewhere jl7| : for 
the charm quark and top quark, electromagnetic spin-spin interactions become important 
and influence t vq and A in such a way as to bring into closer agreement the classical and 
quantum results. 

VII. CONCLUSIONS 

The salient logic in this paper's result is the following. The measurement of a s (M z o) 
and the one-loop A obtained from the so-called 'modified minimal subtraction scheme' of 
renormalization theory [7] predict a s (2mt,c 2 ). From this we may use the classical ballistic 
trajectory lifetime of the VQAP to compute the t vq and action A, obtaining A ~ \% in 
approximate agreement with quantum mechanics. From this we may go further. 

Since % is a more universal and fundamental parameter than a s , % intuitively would seem 
to be the governing parameter in the action equation Eq. ffl5l) . 

If ot s {2m h c 2 ) equalled 3/16 then A would exactly equal Setting Eq. (j3J equal to 3/16, 
Q 2 equal to (2m(,c 2 ) 2 , arid solving for A yields A = 0.106 GeV instead of the standard 0.093 
GeV. With this value of A, Eq. ® gives a s (91.2 GeV = M z o) = 0.121. This is only 2% 
different from the measured value upon which the accuracy of QCD depends, and is within 
the statistical uncertainty in a s (M z o) quoted in Ref. Q. 

We now have a mathematical link between the measured a s (M z o) and the action integral 
for QVAPs from the Uncertainty Principle, \U. This logical sequence is equally valid in 
reverse. We may take as starting point the action A and infer that the action integral is 
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what governs the value of a s (Mzo). This reverse argument from A = hh through Eqs. (fT5|) . 
f fl~2]) . and (j4]) to a s (M z o) is the derivation mentioned in this paper's title. 

The Uncertainty Principle is more general than QCD, and the reason that this derivation 
works is that the fluctuation lifetime At is a different condition than the semiclassical model 
lifetime which is derived from phenomenological QCD. Because these two theoretical life- 
times agree in the special bottom quark case, we can now consider the QCD input parameter 
a s (M z o) to be derivable from our theories. 

This good agreement between the classical trajectory lifetime and the quantum uncer- 
tainty lifetime at the key mass-energy of the bottom quark is surprising, but it may have 
a simple physical explanation: if the vacuum creates these particles in motion at v ~ c, 
then their de Broglie wavelengths A = h/p should be small, and wave packets that are small 
relative to R ma x would represent the particles well. Ballistic dynamics of point masses then 
would serve as a good approximation for the particle motions and the agreement of the 
ballistic and quantum mechanical timescales would be accounted for. The length scale of 

n 

any VAP is usually characterized in standard literature by assuming that v ~ c p. 

The semiclassical model described herein may be used to produce predictions for experi- 
ments involving VQAPs that are interacting with other particles, instead of the unobserved 
VQAPs that were modelled in this paper. 

Acknowledgments 

The author is grateful to NASA's Goddard Space Flight Center for support during this 
research. 



[1] H. Fritzsch, M. Gell-Mann, and H. Leutwyler, Phys. Lett. B 47, 365 (1973). 

[2] D. J. Gross and F. Wilczek, Phys. Rev. Lett. 30, 1343 (1973). 

[3] H. D. Politzer, Phys. Rev. Lett. 30, 1346 (1973). 

[4] W. Lucha, F. F. Schoberl, and D. Gromes, Phys. Reps. 200(4), 127 and see p 150 (1991). 

[5] W. Greiner, Quantum Electrodynamics of Strong Fields (Plenum, New York, 1983). 

[6] J. J. Sakurai, Advanced Quantum Mechanics (Addison- Wesley, Redwood City, Calif., 1967). 

[7] S. Bethke, J. Phys. G: Nucl. Part. Phys. 26, R27 (2000). 

10 



[8] E. Eichen and et al., Phys. Rev. Lett. 34, 369 (1975). 

[9] E. Eichen and et al., Phys. Rev. D 17, 3090 (1978). 
[10] D. B. Lichtenberg, The Standard Model of Elementary Particles (Bibliopolis, Napoli, 1991). 
[11] C. Amsler and et al. (Particle Data Group), Phys. Lett. B667, 1 (2008). 
[12] J. D. Jackson, Classical Electrodynamics, 2nd ed. (Wiley, New York, 1975). 
[13] G. A. Korn and T. M. Korn, Mathematical Handbook for Scientists and Engineers, 2nd ed. 

(McGraw Hill, New York, 1968). 
[14] C. Lanczos, The Variational Principles of Mechanics, J^th ed. (Dover, New York, 1986). 
[15] F. Brau, Phys. Rev. D 62, 014005 (2000). 

[16] B. Sheehy and H. C. von Baeyer, Amer. J. Phys. 49(5), 429 (1981). 
[17] D. Batchelor, in preparation (2010). 



11 



